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Abstract. - We show that the quantum phase transition arising in a standard radiation- 
matter model (Dicke model) belongs to the same universality class as the infinitely-coordinated, 
transverse field XY model. The effective qubit-qubit exchange interaction is shown to be pro- 
portional to the square of the qubit-radiation coupling. A universal finite-size scaling is de- 
rived for the corresponding two-qubit entanglement (concurrence) and a size-consistent effective 
Hamiltonian is proposed for the qubit subsystem. 



Quantum phase transitions (QPTs) are associated with a dramatic change in the phys- 
ical properties of a system at zero temperature when a parameter varies around its critical 
value. It is well-known that very different systems can exhibit similar behavior in this criti- 
cal regime, giving rise to the concept of universality. Enlarging a given universality class by 
the addition of systems from very different areas of physics, is a very important step toward 
unifying our understanding of the basic physics underlying apparently disconnected complex 
phenomena. Recently there have been studies of light-controlled condensed matter systems 
displaying QPTs with atoms in extreme one-dimensional confinements [1] , ions driven by prop- 
erly tuned and pulsed light [2] and fermionic atoms in optical supcrlattices [3] . Fully quantum 
mechanical models of radiation-matter systems are also being considered, and are important 
for several reasons: Scalable and distributed quantum information processing (QIP) devices 
will demand the integration of matter quantum bits (qubits) such as atoms, trapped ions, 
semiconductor quantum dots or SQUIDs with photons. In addition, the capability of photons 
to control and modify the coupling between physically distant qubits makes them appropriate 
for manipulating and transferring quantum information. 
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QPTs themselves are beginning to be studied using concepts from quantum information 
theory. Special attention has been paid to the connections between subsystem entanglement 
and critical phenomena [4-6] . There is an obvious correspondence between condensed matter 
magnetic systems comprising TV spin-1/2 particles, and quantum registers with N qubits. 
Most of the quantum information-based analysis on QPTs performed to date corresponds to 
critical magnetic systems where the exchange coupling between spins is fixed from the outset 
by Nature. However there is also interest in studying QPTs where the spin/qubit couplings 
can be controlled externally, using light for example. One of the most famous radiation-matter 
systems giving rise to a QPT is the Dickc model [7,8]. In this model, N two-level systems 
(qubits) are coupled to a single radiation mode (cavity). Although no direct coupling exists 
between qubits, a superradiant phase arises as the qubit-cavity coupling strength increases. 
Quantum chaos and entanglement in the Dicke model have been studied recently using spin- 
boson transformation techniques [9, 10]. However a full understanding of the entanglement 
(concurrence) finite-size scaling behavior and the associated universality class, are still lacking. 

Here we show that the reduced dynamics of a qubit system embedded in a photon field, as 
described by the Dicke model, displays magnetic-like behavior and corresponds to the same 
universality class of quantum critical phenomena as an infinitely-coordinated ferromagnetic 
system. The qubit-cavity detuning determines an effective magnetic field while the square of 
the qubit-cavity coupling strength yields the effective qubit-qubit 'exchange' interaction. As a 
result, we are able to establish an unambiguous quantitative connection between the associated 
scaling exponents and effective qubit Hamiltonian in two quite different physical systems. 
Although Lambert et al. [10] have already noted similarities between the Dicke and collective 
XY model in terms of scaling exponents for the concurrence, the present work establishes a 
systematic link between these models and lays these similarities on a firm theoretical footing. 

We are interested in a system formed by N identical two-level atoms (qubits) coupled to 
a single-mode cavity, described by the Hamiltonian 

N \ ^ 

where a^ (a) is the creation (annihilation) operator for single-mode cavity photons, ai^z = 
\li){U\ - \0^){0^\, a] = |1,)(0,|, a, - |Oi)(li| with |0,) and |1,) {i = 1,2, ...,7V) being the 
ground and excited states of the z'th qubit, respectively, and A is the qubit-cavity coupling. 
Note that we do not make a rotating wave approximation. We take the energy of a confined 
photon as the unit of energy. Resonance between qubit energy levels and photons corresponds 
to e = 1. This Hamiltonian can also be written as 

2A 

H ^a)a + eJ^ + ^={o) + a)J^ (2) 

ViV 



where Jz = o Si=i ^j,^ *^^ Jx — h'^i=ii"i + '''i)- We denote this latter Hamiltonian as 
H = Ho + V where Hq = a'^a and V ^ eJ^ + -%(a^ + ajJ^. 

V-/V 

We shall be concerned with physical qubit quantities such as 'magnetization' and en- 
tanglement. These quantities can be represented by a generic qubit operator Q. Hence in 
thermodynamic equilibrium at inverse temperature (3 — {kgT)^^ any qubit expectation value 
can be expressed as 

Tr,,[QTrp4e-n] ^ Tr,,[Qe~0"^^] 
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where Tr^f, {Tvph) denotes a trace over the qubit (cavity) subsystem and the last equahty is 
the formal definition of an effective qubit Hamiltonian Hqi,. In order to eliminate the cavity 
degrees of freedom and obtain a size-consistent Hqi, (i.e. a reduced qubit Hamiltonian which 
scales with the system's size) we follow the procedure sketched out by Polatsek and Becker [11]. 
The resulting effective qubit Hamiltonian can be written as 



H. 



qb 



il^(^^-PiLo + V)) 



13 



ph 



(4) 



where the thermal averaging is carried out with respect to the cavity degrees of freedom, 
i.e. {0)ph = Trph{Oe~P°- °')/Trph{e~^°' "); Lq is the Liouville operator associated with iJo 
and defined by LqO — [a^a, O] for any operator O. A cumulant expansion is well-suited to 
calculating a size-consistent TJ^b, yielding 



n—l 



(5) 



where Hn represents the qubit-operator corresponding to the n'th cumulant. Cumulants 
coming from eq. |0J are themselves operators which act on qubit states, hence any expression 
involving cumulants must be properly symmetrized. The lowest-order cumulants are: 



Hi = eJ, 



H2 = 



2\ 



{2h{(3) + l)Jl 



(6) 



where h{0) = [e^ — 1)^^ is the Bose factor which determines the average photon number in 
an isolated cavity at temperature T . Due to the superradiance phenomenon which occurs in 
the Dicke model, the actual photon number has very different behavior on either side of the 
quantum and/or thermal phase transition. It can be seen that Hi corresponds to the non- 
interacting qubit system while H2 gives rise to the lowest order effective qubit-qubit coupling. 
Every cumulant contains a term of the form 



{VL-^V)pH = 



2A 



[{i + {^ir)h{f3) + i]ji 



(7) 



which enables us to perform a closed- form summation in eq. 10 over an infinite set of contribut- 
ing terms. Up to second order in V , this renormalization summation yields a temperature- 
dependent effective Hamiltonian given by 



H'^^{fl)^eJ, 



2A 



[1 



I3{h{(3) + 1) 



— 1 /2 



(8) 



At zero temperature (i.e. (3 — > 00) this effective Hamiltonian is a special case of that 
studied in ref. [6] corresponding to a ferromagnetic infinitely-coupled XY model in a transverse 
magnetic field with a QPT at 2 Ac = 1 in the thermodynamic limit (i.e. N — s- 00). Our 
analysis therefore provides a novel demonstration of the well-known QPT arising in the Dicke 
model at Ac ~ 1/2 [8]. The broken symmetry associated with the phase transition is related to 
parity for the Dicke model as given by e"*^" a+j^+^) ^ while for the effective H, model this 
becomes e"('^'"'"~\ In contrast to a previously employed continuous variable representation 
[9, 10], this effective anisotropic XY model Hamiltonian H\' provides an ideal starting point 
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Table I - Critical exponent b at zero temperature for the finite-N scaling of the x-magnetization at 
Ac (first row), z-magnetization at Ac (second row), difference between the concurrences at Ac for the 
thermodynamic limit and for N qubits (third row), difference between the maxima of the concurrences 
for the thermodynamic limit and for N qubits (fourth row), and the difference between the coupling 
strength for N qubits at maximum concurrence and Ac (fifth row). N < 35 for the Dicke model and 
100 < N < 500 for the anisotropic XY model. 



f{N) ~ iV" 


Dicke 


XY 


N 


-0.35 ±0.01 


-0.34 ±0.01 


1 
2 


N 




-0.54 ±0.01 


-0.55 ±0.01 


Coo (Ac) "Civ (Ac) 


-0.26 ±0.01 


-0.30 ±0.01 




-0.28 ±0.03 


-0.30 ±0.03 


Ajv — Ac 


-0.65 ± 0.03 


-0.66 ±0.03 



for studying the QPT of the Dicke model since we can follow the qubit subsystem's quantum- 
mechanical evolution as a function of the effective qubit-qubit coupling. Furthermore, this 
mapping onto the infinitely-coordinated XY model which results from our approach, yields a 
deep understanding of the scaling properties for finite N as well as the mean-field behavior of 
the Dicke model in the thermodynamic limit. Fig. ^ shows the ground-state energy (at zero 



o 




Fig. 1 



Fig. 2 



Fig. 1 - Ground state energy and z-magnetization (inset) as a function of A for A'' = 15 qubits, as 
obtained from the full Dicke Hamiltonian (x) and from the effective Hamiltonian (o). 

Fig. 2 - Non-rescaled concurrence as a function of A for the XY model and for the full Dicke Hamil- 
tonian (inset) for Y=15 and 35. 



temperature) of a Dicke system with N qubits, as a function of the qubit-radiation coupling A. 
For simplicity, we only consider the case of resonance from now on, i.e. e = 1. For comparison, 
the ground-state energy produced by the effective Hamiltonian in eq. ^ is also included in 
the plot. The agreement between both ground-state energies is excellent even for low N . The 
qubit subsystem's 'z-magnetization', which is the order parameter for this QPT, undergoes a 
quantum phase transition at Ac = 1/2 as seen in fig. Q] (see inset) for both the exact Dicke 
model and the effective Hamiltonian. This indicates that the physics of the QPT in the Dicke 
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Fig. 3 - Dependence of Xff — \c (lower curves) and C^ — Cfi (upper curves) on A'' for H (A) and 
H-^' (x). A^ is the value of A for which the concurrence is maximum, while Cff is the maximum of 
the rescaled concurrence. 

Fig. 4 - Finite-size scaling for the rescaled concurrence in the Dicke model. Points corresponding to 



different A'^ collapse on the same curve. Inset: Finite-size scaling for the x-magnetization rrix 



model is indeed captured by the effective Hamiltonian of eq. ((HJ . 

As a more stringent test of the effective Hamiltonian in eq. IjSJl, we now compare other 
thermodynamic as well as quantum features. An important ground-state quantum feature 
is the entanglement between two qubits. We adopt the concurrence Cnr as a measure of 
this entanglement [12]. We also study the rescaled concurrence C = NCmr m order to 
obtain non- trivial information in the thermodynamic limit where Cmr — Q. In contrast 
with thermodynamic results that only involve the calculation of the partition function and 
averages obtained from it, the concurrence probes the internal structure of the ground-state 
wavefunction in a more detailed way. In fig.|21we display the evolution of the concurrence as 
a function of A, for both the exact Dickc model and the effective Hamiltonian of eq. |(SJ). The 
qualitative form of the concurrence for both models is very similar, even though the effective 
Hamiltonian overestimates the maximum value. This overestimation is easily understood: 
the effective Hamiltonian drops higher-order interaction terms which incorporate additional 
interactions among the qubits and would hence modify the ground-state wavefunction. In the 
thermodynamic limit, the rescaled concurrence displays a singularity at A = Ac in agreement 
with the findings in refs. [6, 10]. However, at Ac the rescaled concurrence goes to 1 for the XY 



model while it goes to 1 



V2 



for the Dicke model. 



As a final test of the reduced Hamiltonian of eq. , we consider the finite-size scaling of 
the 'magnetization' as well as the concurrence. Results can be compared with those reported 
in ref. [13] and show excellent agreement between the predictions of the Dicke model and the 
effective Hamiltonian. 

At zero temperature, the energy and the z-magnetization are identical for both systems 
in the thermodynamic limit. Physical quantities depending on the ground state structure, 
such as the x-magnetization and concurrence, display some differences. However, the scaling 
exponents agree for both models as can be seen in fig. O and table Q] Even for small N in the 
Dicke model, the quantities X^ — Ac and C^ — C^^ can be seen to lie on well-defined straight 
lines. For the XY model, large N results lie on straight lines which are parallel to the Dicke 
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model results. The Dicke model results approach the thermodynamic limit much faster than 
those for the XY model. This is directly related to the fact that terms have been thrown 
away in eq. ((Sj) which include interactions between three and more spins, thereby enhancing 
the mean-field nature of this system. 

The excellent agreement between the critical exponents for the two models provides clear 
evidence that they both belong to the same universality class. It is well known that in both 
Hamiltonians, each qubit is coupled to every other qubit with the same strength and hence 
no typical lengthscale can be associated with them. Botet et al. [13] proposed relating the 
infinitely-coordinated magnetic system to short-range interaction models, and hence were 
able to perform a systematic study of finite-size scaling functions for the x-magnetization and 
gap energy in the case of a Hamiltonian such as H\' . Following this line of reasoning, we 
show in fig. 21 that this scaling hypothesis is valid for the Dicke model as well. In particular 
there exist scaling functions in the critical region |A — Ad — > 0: the concurrence behaves 
as Coo(Ac) - Cn[X) = |A - Ac|""^Fc(Ar|A - A^'') with F^ a function of N\X - A^r only. 
When a; — > CX3, Fc[x) -^ constant and Coo(Ac) — Coo(A) ^ |A — Ac|°™^ while for x ^ 0, 
Fc{x) -^ x-'^-^f^" and Coo(Ac) - CAr(Ac) ~ N^'^-^f/". In ref. [14] the concurrence in the XY 
model {N — > cx), eq. (jHJ) has been expHcitly calculated yielding 



l-A/l-lf)' , 0<A<Ac 



C(A) = V^ A^ (9) 

l-\/l-(#) , Ae<A<cx) 



with Ac = -^. This implies that, close to Ac, 1— C^ |A — Ad^. If we assume for the moment 
that both models, Dicke and infinitely coordinated XY, belong to the same universality class, 
we can then borrow the result amf = 5 from the XY concurrence and apply it to the Dicke 
concurrence. This exponent has also been found for the full Dicke model in ref. [10] using 
a spin-boson transformation scheme. The exponent — ^ has been numerically obtained (see 
table 0} yielding 1/ = 1.88. The collapse of all the data on a single curve for each region 
(coupling below or above the critical coupling strength) is a verification of the existence of 
scaling behavior of the concurrence in the Dicke model. In fig. 21 (inset) the x-magnetization 
multiplied by ]A — Ac]"^/^, is plotted as a function of A^] A — Ac] ^■'*'^. The critical exponent in 
the argument of the scaling function is 1.43, which is only slightly smaller than that found for 
the same exponent in the XY model case (1.5). This similarity in scaling exponents reinforces 
the main message of this paper, that the scaling behaviour of the Dicke model is intimately 
related to that of the anisotropic infinitely-coordinated XY model. 

In summary, we have shown that a radiation-matter (i.e. Dicke) model can be mapped on 
to an infinitely-coordinated spin/qubit system with an effective magnetic field determined by 
the detuning between qubit and cavity, and a long-range exchange interaction given by the 
square of the qubit-cavity coupling strength. An immediate and interesting consequence is that 
the exchange interaction can be controlled experimentally simply by changing the system's 
temperature or equivalently the mean number of photons in the cavity. We have also shown 
that the Dicke model pertains to the same universality class as other infinitely-coordinated 
systems, such as the anisotropic XY model in a transverse field. We have demonstrated that 
the finite-size scaling hypothesis works well for the Dicke model, and we therefore postulate 
that it is also likely to work well for a range of generalized Dicke models corresponding to a 
diverse range of radiation-matter experimental setups. Temperature-dependent phenomena 
can also be studied using our size-consistent, effective qubit-qubit system approach. 
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